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Swift captures the spectrally evolving prompt emission of 
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ABSTRACT 

The origins of Gamma-ray Burst prompt emission are currently not well understood 
and in this context long, well-observed events are particularly important to study. 
We present the case of GRB 070616, analysing the exceptionally long-duration multi- 
peaked prompt emission, and later afterglow, captured by all the instruments on-board 
Swift and by Suzaku WAM. The high energy light curve remained generally flat for 
several hundred seconds before going into a steep decline. Spectral evolution from hard 
to soft is clearly taking place throughout the prompt emission, beginning at 285 s after 
the trigger and extending to 1200 s. We track the movement of the spectral peak en- 
ergy, whilst observing a softening of the low energy spectral slope. The steep decline in 
flux may be caused by a combination of this strong spectral evolution and the curva- 
ture effect. We investigate origins for the spectral evolution, ruling out a superposition 
of two power laws and considering instead an additional component dominant during 
the late prompt emission. We also discuss origins for the early optical emission and 
the physics of the afterglow. The case of GRB 070616 clearly demonstrates that both 
broadband coverage and good time resolution are crucial to pin down the origins of 
the complex prompt emission in GRBs. 

Key words: gamma-rays:bursts 



1 INTRODUCTION 

The prompt emission mechanism for Gamma-ray Bursts 
(GRBs) is commonly attributed to internal shocks due 
to collisions of shells of different Lorentz factors ejected 



* This paper is dedicated to the memory of Dr. Francesca Tam- 
burelli who died during its production. Francesca played a fun- 
damental role within the team which is in charge of the develop- 
ment of the Swift XRT data analysis software at the Italian Space 
Agency's Science Data Centre in Frascati. She is sadly missed. 



from the vicinity of a compact object (jRees fc Meszaroj 
1 1994 ). Later afterglow emission comes from an external 
shock as the GRB blastwave dece lerates through interac - 
tion with the surrounding medium IjMeszaros fc Reeslll997l ). 
Though the general picture appears applicable to most 
GRBs, the details are far from understood and other mod- 
els, in part icular magnetised flows, have also bee n pro- 
posed (e.g. Meszaros. Rees fc Papathanassioul Il994l : lUsovl 
ll994l : lKumar et al.ll2007l : IZhang||2007l v 



The advent of the Swift mission (JGehrels et al.l |2004 ) 
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has revealed additional f eatures, for ex a mple s teep de- 
cays and X-ray flares ([Burrows et alj l2005al ). whose 
properties are consis t ent with an intern a l origin (e.g. 
Tagliaferri et all l2005l: iNousek et"aLl l2006l: lO'Brien et al] 



20061 : IChincarini et al.1 l2007l : iFalcone et alJ 120071 ). The 
steeply decaying phases that directly follow both prompt 
emission and X-ray flare s are usually interpreted as due to 
the curvature effect (e.g. iKumar fc Panaitescull2000l ) where 
high latitude emission is delayed with respect to that on- 
axis. However, significant spectral evolution is not expected 
in this model, and new mechanisms must be invoked to ex- 
plain those observations that do show spectral evolution dur- 
ing the steep decay phase (Zhang, Liang & Zhang 2007). 

Of particular interest are the longest duration GRBs 
in which the relationship between various possible early 
emission components can be studied. Observationally, very 
few GRBs are detected in 7-rays for more than 400 sec- 
onds as quantified using the T90 parameter. For example, 
appro ximately 0.5% of the BATSE sample meet this cri- 
teria IjPaciesas et al.l Il999r ) . Some of the very long dura- 
tion GRBs in clude those with a pre-cursor or a late flare 
fe.g. iPrice eT al. 2002; Ni castro et aLlbOO^ ). Interestingly, 
the very long GRBs also include several FRED-like (fast 
rise, exponential decay) events which can be sp ectrally quite 
soft l|Giblin et al.ll2002i : liiT't Zand et al.ll2000l ). 

In the Swift era it remains true that very few GRBs 
have a T90 > 400 s. The longest is GRB 060218, which 
is an unusually l ong, lo w luminosity, spectrally soft GRB 



IIGampana et al. 2006b). In the first Swift BAT cata- 



logue (|Sakamoto et al 



20071 ) there are only three others: 
GRB 060929, which has a spectrally soft, late giant flare 
jPalmer et al. 2006): GRB 070129. for which BAT triggered 



on a precursor (|Krimm et aLll2007l ): and GRB 070616. 

Here we present data for the case of GRB 070616, in 
which the prompt emission shows a very complex multi- 
peaked structure, leading to one of the longest prompt emis- 
sion durations ever recorded. We take advantage of extensive 
coverage of s uch a long burst by the Swift Burst Alert Tele- 
scope (BAT; 'Barthelmv ct all boosh . the X-Ray Telescope 
(XRT;Purrows ct al. 200510^ and the UV/Optical Telescope 
fUVOTT TRoming et al.ll2005l). C ombining data from Swift 
and Suzaku ( Mitsuda et al.ll2007r ) we study the evolution of 
the prompt emission spectrum, following the temporal vari- 
ability of the peak energy, and study the development of the 
afterglow component. 

In Section 2 we describe the Swift and Suzaku obser- 
vations. In Section 3 we outline the X- and 7-ray temporal 
characteristics. In Section 4 we model the X-ray to 7-ray 
spectra and investigate spectral evolution. We discuss the 
discovery of an optical transient in the UVOT images in 
Section 5 and go on to model the spectral energy distribu- 
tions in Section 6. In Section 7 we present possible interpre- 
tations for our findings, and discuss them in the context of 
the blastwave model and in comparison with other GRBs, 
with summary and conclusions in Section 8. 



2 OBSERVATIONS 

On 2007 June 16 at 16:29:33 UT (hereafter To), the Swift 
BAT triggered on and located GRB 070616 (trigger=282445. 
Starling et al. 2007a). This was an image trigger, where the 



GRB is detected by forming an image from the collected 
counts which is then searched for a new point source. Swift 
slewed immediately to the BAT location. T90 (15-350 keV) 
is 402 ± 10 s (estimated error including systematics), which 
is among the longest T90 values for BAT GRBs. 

The XRT began observing the field at 16:31:44 UT, 131 
seconds after the BAT trigger; a bright, fading and uncata- 
logued X-ray source was clearly detected. XRT observed the 
source for a total exposure of 144.5 ks out to To-|-3.7xlO^ 
s, beginning in Windowed Timing (WT) mode at To-|-137 s 
and continuing in Photon Counting (PC) mode at To -1-976 s. 
Using 899 s of overlapping XRT PC mode data and UVOT 
V band data, we obtain an astrometric ally corrected X-ray 
position (using the USNO-Bl catalogue. ICoad et aLll2007bh 
of RA (J2000) = 02h 08m 36.59s; Dec (J2000) = -h56° 56' 
43.8", with an error of radius 2.4 arcseconds (90% contain- 
ment). The error circle contains the bright star USNO-Bl.O 
1469-0076513. 

The UVOT took a finding chart exposure of 100 seconds 
with the White filter starting at To -1-142 s, and continued 
with a long V band exposure followed by cycling through 
all seven filters. The count rates in V and U band exposures 
between 250 and 1050 seconds after the trigger show a sig- 
nificant excess (at ^3cr) with respect to late observations 
(>1 day after the trigger) at the position of the USNO-Bl.O 
star. We discuss this further in Section 5. 

The Suzaku W ide-band All-sky Monitor (WAM; 
lYamaoka et al.l 20051). one element of the Hard X-ray De- 
tector (HXD: lTakahashi et al.ll2007^ with total bandpass 50 
keV to 5 MeV, detected GRB 070616 at 16:31:50 UT, 137 
s after the Swift BAT trigger. The observations cover 140 s 
of the GRB emission, truncated by passage into the South 
Atlantic Anomaly. With these data alone we put a limit on 
the GRB duration of T90 of > 112 s over 150-700 keV, con- 
sistent with the BAT (15-350 keV) T90. 

The Swift data were processed with the standard Swift 
data reduction pipelines, and spectra and light curves were 
extracted with xselect. All XRT spectra are grouped such 
that a minimum of 20 counts lie in each bin, and X-ray 
light curves have a minimum of 15 source region counts per 
bin. We used a source extraction region of box width 40 
pixels for XRT spectra in WT mode (1 pixel = 2.36"). In 
the first orbit of PC mode data, emission from the source was 
piled up, and spectra were extracted using a 20 pixel radius 
annulus excluding the inner 4 pixels in radius; a correction 
for this was made to the ancilliary response file. Thereafter 
a circular source region of radius 17 pixels was employed. 
Background spectra were extracted from a nearby source- 
free region of r adius 86 pixels. Spectral fitting was performed 
in Xspec (JArnaud 1996 ) using version 008 response files. 
Errors are given at the 90% confidence level throughout, 
unless otherwise stated. 



3 TEMPORAL CHARACTERISTICS OF THE 
PROMPT EMISSION 

Light curves for the prompt emission in all observed bands 
are shown in Fig. [l] The 7-ray light curves show no strongly 
peaked or variable emission around Tq, as is common for 
BAT image triggers such as this. BAT triggered instead on 
a gradual rise which lasted approximately 100 s before the 
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Figure 1. Light curves for the Suzaku WAM and all Swift instru- 
ments are shown in count s~^. The instrument and, where appro- 
priate, the energy range used in keV, are given in the y-axis labels. 
We show here the WAM 1 s light curve and the mask-weighted 
1 s light curves for the BAT in 4 individual energy ranges. The 
UVOT observed in all seven of its filters, but we only show the 
V band here (500-560 nm), in which the majority of the detec- 
tions were obtained, for clarity; USNO-Bl.O 1469-0076513 has 
been subtracted. 



first and strongest peak, centred at To-|-120 s. Thereafter 
multiple blended peaks continue the prompt emission out to 
To-|-500-600 s. At tliis point the 7-ray emission appears to 
return to the count rate at whicli it began at To, wliilst tlie 
X-ray emission begins a steep decay lasting until To-|-1200 s 
(Fig. (2)1 , and the V band optical emission continues a steady 
rise (Fig. [TJ. 

The X-ray observations cover much of the time period 
over which BAT could detect the source, and we constructed 
a joint BAT-XRT light curve to compare emission from the 
two bands. This was done by performing joint power law 
fits to fine time-sliced BAT-XRT (WT) spectra, and ex- 
trapolating the BAT spectra to the XRT energy band. The 
main peaks in the multipeaked prompt emission are found 
to be temporally coincident, strongly suggesting that the X- 
ray and 7-ray emission come from the same component: the 
GRB prompt emission. In contrast, the V band optical data 
show a slowly rising light curve from To -1-250-1000 s which 
is not mimicked in any of the high ener gy bands (Fig, [l]). 

We performed a lag analysis (e.g. iNorris et al.l Il996l ) 
over two intervals in the BAT light curve: ~70 s covering 
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Figure 2. The XRT 0.3-10 keV count rate light curve, plotted in 
log space. 



the largest peak, and ~155 s covering the flat multipeaked 

structure following the largest peak. We use the four BAT 

channels 15-25 keV (channel 1), 25-50 keV (2), 50-100 keV 

(3) and 100-350 keV (4) and 256 ms time binning. For the 

first interval we find the following lags: 

3^1 = 220t?5o ms ; 4^2 = 410±200 ms 

and in the second interval the lags have decreased to: 



f-80 
-112 



The time resolution of the analysis prevents accurate mea- 
surements of lags <256 ms, causing the relatively large errors 
reported here. 

We attempted to model the underlying temporal be- 
haviour of the prompt X-ray emission by fitting a multiply 
broken power law model to the XRT light curve, as is often 
done for complex GRB light curves. We find the following 
temporal slopes, q (F oc t~"), and break times Tbk: Qi = 



-(0.08: 



■0.01. 

18 

'24 



up to Tbk.i = 524 



+3, a2 



4.87 



-1-0.11 



up to Tbk,2 



999:1:24 ^id Q^3 > 10 up to 1200 s (the steepest segment of 
the decay has a slope so large when fit with a simple power 
law it is difficult to accurately measure). Between the steep 
slope and the first probable afterglow emission. To -I- 1130 to 
To-|-5800 s, we measure Opiateau = 0.5lo 2 using a power law 
fit to the last point on the steep decay and the second orbit 
of PC mode data. The afterglow temporal slope from >4600 
s using this method is 04 = 1.46±0.05. 

For the XRT and BAT we obtained hardness ratios as 
a function of time (Fig. [3]). The BAT hardness ratio (15- 
50 keV/50-100 keV) remained approximately constant until 
To -f 285 s when the spectra softened significantly over the 
remainder of the 7-ray observations. The XRT hardness ra- 
tio (1-10 keV/0.3-1 keV) shows the same behaviour, with 
the spectral evolution beginning at approximately To-|-500 
s (see also Fig[lO]in Section 7.1). In both bands the emission 
begins with a hard spectrum, evolving to become softer. The 
source is spectrally harder at higher count rates, as shown 
in Fig. |4] In the lower panel of Fig. |4] we show that flares 
themselves are not driving this correlation, although they 
are consistent with it, but that there is an overall trend in 
the hardness ratios with only the smaller scale changes cor- 
responding to hardening within individual flares (individual 
flares can be seen as loops, following the colour scale from 
red through to blue which represents the passage of time). 
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Figure 3. Count rate (CR, upper panel) and hardness ratio (HR, 
lower panel) against time for XRT (1-10 keV/0.3-1 keV, 10 s bins, 
grey) and BAT (15-50 keV/50-100 keV, 10, 20 or 50 s bins, colour 
scale) . 



4 SPECTRAL FITTING 

4.1 X-ray 

To investigate the spectral evolution seen in Figs. [3] and HI 
we began by characterising the XRT X-ray spectrum for the 
prompt emission interval. All results are given in Table [T] To 
test for the presence of intrinsic absorption, we extracted a 
time-averaged spectrum over the time range To -I- 137-963 s. 
We fit the spectrum with a power law model with absorption 
set initially to the Galactic value of A^n .Oai = 0.35 x 10^^ 
cm~^ (LAB survey. iKalberla et al.ll200a) and then allowed 
to go free. The inclusion of a small column of intrinsic ab- 
sorption, amounting to 5 x 10^" cm~^ at 2: = 0, is required 
according to the F-test (F-statistic = 6xl0~^^, Table[T}. 

We then extracted X-ray spectra in 100 s time bins 
from To-l-137 to To-|-837 s, and two further longer time 
bins of sizes 126 s and 224 s to extend coverage to the 
end of the steeply decaying phase (To-|-1200 s). We fit these 
spectra with a power law plus Galactic-|-intrinsic absorption 



(amounting to a total column of A^h = 0.4 x 10^^ cm~^). We 
also allowed the intrinsic absorption column to vary to test 
for any temporal changes which may be caused by e.g. ioni- 
sation by the GRB emission, but found no variation confirm- 
ing that the time-averaged A'h, intrinsic that we have applied 
is appropriate throughout. We observed the power law slope 
to change from F ~ 1.1 to F ~ 2.3 in ~800 s, indicating 
that a spectral softening occurred in the X-ray band. This 
begins after the To -1-337-437 s segment, and corresponds to 
the time of the change in the X-ray hardness ratio (see Fig. 

EJ. 

We also fit a broken power law model to the same data, 
initially with Galactic absorption only. The broken power 
law model provides statistically equal or better fits than 
the absorbed power law model in most cases, but this may 
be a symptom of the need for intrinsic absorption, and/or 
possibly additional spectral curvature. Adding 5 x 10^° cm~^ 
of intrinsic X-ray absorption does not significantly improve 
the majority of the broken power law fits. Although it is 
difficult to statistically distinguish between a power law plus 
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Figure 4. Hardness ratio against count rate for XRT (upper) 
and BAT (lower). The evolution of both quantities with time is 
shown by the colour gradient in the BAT panel, going from red to 
blue as time increases from To to To -1-650 s. This colour gradient 
is identical to that shown in Fig. \3\ The general trend of harder 
spectra at higher count rates is visible, and within this we see 
individual loops of soft-hard-soft evolution. Each loop shows the 
effect of the larger individual flares. 



intrinsic extinction and a broken power law, we find that 
the break energies and power law slopes in the latter model 
move seemingly randomly in time whereas the steepening 
in the single power law model matches the hardness ratio 
behaviour very well. In general, we expect that there will 
be some material within the host galaxy that lies in front 
of the GRB and alon g our line-of-sight caus ing detectable 
X-ray absorption (e.g. lCampana et al.ll2006al ). We therefore 
do not consider the broken power law model further in fits 
to the XRT data. 

The X-ray afterglow during the second orbit of PC 
mode data (To -1-4608-6976 s) shows a soft spectrum which 
is well fitted by a single absorbed (Galactic -|-intrinsic as 
above) power law with F = 2.53J:E;;i^ (xVdof = 11/13). No 
spectral variation is seen through the afterglow phase when 
comparing the above time range with data from To-l-10 ks 
onwards which are well fit with a power law of photon index 



2.62^ 



(xVdof =14/16). 
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Table 1. X-ray spectral fits, using the XRT energy range 0.3—10 
keV, and over varying time ranges covering the prompt phase: 
To-|-137-1200 s. We test the single power law model with various 
absorption columns applied. Fstat refers to the F-test statistic 
when comparing the fits with absorbing column allowed to vary 
against those fixed at the Galactic column of 0.35x10^^ cm~^. 
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0.35 
0.35 
0.35 
0.35 
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0.35 
0.35 
0.35 
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gq+0.02 

-, iq+0.02 

^•^•^-0.03 
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4.2 X-ray to 7-ray 

We now extend the above analysis to include the higher en- 
ergy BAT spectrum. Firstly we extracted a time-averaged 
BAT-XRT spectrum covering the prompt emission phase 
from the start of XRT observations to the end of source visi- 
bility with the BAT (137-963 s). We fit this with three mod- 
els: a power law, broken power law and the Band function 
l|Band et al.lll993h . The normalisation for each instrument 
is always tied such that they are equal. In the Band function 
fits we fixed the high energy po wer law slope to /? = 2.36 
as found in the Suzaku spectrum IjMorigami et al.ll2007l . and 
adopting Fi, oc i^~ ), since the energy range with slope /3 
is expected, at least initially, to lie above the BAT band. 
Each model was also absorbed by either the Galactic col- 
umn alone, or the previously determined Galactic -|-intrinsic 
column. Again inclusion of the small amount of intrinsic ab- 
sorption improved the power law and Band fits, but was 
not a significant improvement for the broken power law fits. 
Results are given in Table 2. A single power law is a poor 



2,0 












: 


1.8 


— 










- 




T 


: 


1.6 










1 


: 


1.4 


— 






T 




— 


1 2 
















_ 


















- 
















150 


- 










- 


"Z 100 


7 


J 




1 


w'' 50 


^ 




T 






- 


1 










I 


T 


_ 













— 1 — 


- 



100 



200 



700 



300 400 500 600 

time since BAT trigger (s) 

Figure 5. The Band model fitted to the BAT-XRT spectra in 
seven 100 s intervals, where the fitted parameters a (upper panel), 
and Ef) combine to give Spj^ (lower panel) .90% errors on Spj; have 
been calculated using average symmetrical error bars for a and 
Eq. 



representation of the time-averaged joint spectrum, and we 
find the Band function or a broken power law provides a 
better fit. 

To follow the spectral evolution, we time-sliced the BAT 
data into six 100 s bins covering To -I- 137-737 s, correspond- 
ing to some of those extracted for the XRT. We fit the three 
models listed above to the joint BAT-XRT spectra for each 
time bin, this time absorbed by the combination of Galactic 
and intrinsic absorption. The single power law model shows 
the same softening of the spectral slope with time as was 
seen for the time-sliced XRT data. The broken power law or 
the Band function again provide a better fit than a single 
power law model. The broken power law model fits suggest 
a fairly constant break in the middle of the X-ray band at 
^4-5 keV (we note that the fits are biased towards the X- 
ray band where there are the greatest number of counts). 
In the Band function fits the peak energy is derived from 
the free parameters a (low energy power law index) and Eo 
(characteristic energy) and the fixed parameter /3 = 2.36 
(high energy power law slope) using 



£pk = £o(2-Q) 



(1) 



(see Band et al. 1993 for a detailed description of the Band 
function; errors on £pk are calculated using average symmet- 
rical error bars). _Epk can be well constrained and is observed 
to move to lower energies with time from 135 keV down to 
4 keV in ~600 s (Fig. (5] Table 2) , while the spectral slope 
a also varies gradually, softening with time. 

To confirm the validity of the high energy spectral slope 
we fixed in the previous Band function fits and the position 
of the peak energy derived from fits to the Swift data alone, 
we performed power law fits to the Suzaku WAM spectrum. 
The constraints we derive on the power law slopes are con- 
sistent with the value /3 — 2.36 used in previous analysis 
(see Table 2). We then performed joint fits for WAM, BAT 
and XRT. The results allow us to confirm the previous find- 
ings using a broader-band spectrum. In particular, previous 
BAT-XRT fits suggest the peak energy of the spectrum lies 
in the Suzaku energy band at these times (^250 s) and is 
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Figure 6. Spectral energy distribution of the WAM, BAT and 
XRT spectra at To+133-159 s, with best fitting absorbed Band 
function. The peak of the spectrum can be seen close to the 
boundary between the BAT and WAM bandpasses. 



therefore more accurately measured in the WAM-BAT-XRT 
joint fit. Results are reported in Table 2. The absorbed Band 
function fit to the broadband data is shown in the E^Fe plot 
(energy equivalent of uFv) in Fig. [6] 



5 IDENTIFYING THE OPTICAL 
AFTERGLOW 

An optical transient (OT) was not immediately identified in 
the Swift UVOT images owing to superposition with the 
bright {V = 14.4) star USNO-Bl.O 1469-0076513, whose 
colours are consistent with those of a hot star. The blended 
image of the star plus the OT is extended compared to other 
stars in the field. In both the UVOT V and U bands the OT 
was separated from the nearby star by subtracting a prop- 
erly scaled later image from earlier images. For the V band 
subtraction we used the sum of the first and third exposures 
(To-l-250-650 s and To-|-986-1152 s) and a sum of 2 later 
exposures covering To-|-178. 2-185.1 ks. The U band sub- 
traction used a sum of the first 2 exposures (To +704-723 s 
and To -1-866-874 s) and a sum of 3 later exposures covering 
To-|-97. 2-98.1 ks. The scaling was based on the total expo- 
sure times and all images were mod-8 and aspect-corrected, 
and rebinned to 2x2 UVOT sub-pixels (l"xl"). The sub- 
tracted images, of which the V band is shown in Fig. [T] 
reveal the sources which have varied between the considered 
epochs. Due to minor changes in the average point spread 
function (PSF) of the stellar source during the orbit, count 
rate noise in the original images, and the combined effect 
of uncompensated coincidence-loss where the OT and stel- 
lar PSFs overlap, the stellar image could not completely be 
removed and is present as a low intensity peak and ring 
of opposite sign. Some residual of other bright stars can 
thus also be found in the subtraction. The net images in 
the V and U filters differ slightly due to different PSF and 
coincidence-loss effects in each band. These complications 
have only a minor effect on the determination of the OT po- 
sition, which is found from a combination of the subtracted 
images of both bands. We find a best position for the af- 




■^ 




Figure 7. UVOT V band image subtraction revealing the fading 
optical counterpart (negative areas are white, while positive areas 
are shades of grey). The white circle is the OT position and the 
white 'x' marks the centre of the star USNO-Bl.O 1469-0076513, 
as measured from both V and U bands. Contours show both the 
residuals of the nearby star and the OT. The black 'x' indicates 
the position of nearby field star USNO-Bl.O 1469.0076499. 1 pixel 
= 1". 



terglow at RA (J2000) = 02h 08m 36.37s; Dec (J2000) = 
56° 56' 44.1", with an uncertainty of 1.2" (radius, 90% con- 
fidence). Using a 4" radius extraction region and a nearby 
source-free background region, we determine the count rates 
and fiux densities for the OT after subtracting those of the 
star (Table El). The OT is detected in the U and V bands 
at >3a, and at lower significance in the B band. It is not 
detected in the UV filters and has disappeared below detec- 
tion limits in exposures taken after 1147 s in V, 739 s in 
B and 870 s in f/. We do not use the white filter here due 
to high coincidence losses. Detection of this source in the U 
band indicates a redshift less than ~3. 

A limit on the near-infrared magnitude was found us- 
ing the 1.34m ground-based TLS Tautenburg telesc ope us- 
ing ISIS image subtraction (|Alard fc LuptorJ 119981 ) of two 
epochs, finding a conservative limit on the magnitude of 
/ > 19 at 0.3 days after the trigger ( K. Wiersema fc 
D.A. Kann, private communication; see also lKann fc Wilsonl 
|2003). 

To support the argument that the OT is the optical 
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Table 3. UVOT background subtracted measurements and lim- 
its for the optical afterglow candidate, performed with a 4 arcsec 
aperture and with the contribution from the USNO-Bl.O star re- 
moved. All errors ar e la and upper limits 3(7. Galactic e xtinction 
of E{B - y) ~ 0.4 llSchlegel. Finkbeiner fc Davis|[T993) has not 
been accounted for. 



Band 



Tc 



Count/s Flux (/iJy) 



317.8 


135 


0.45±0.42 


113±105 


416.7 


62.6 


2.09±0.61 


530±154 


478.9 


61.5 


2.49±0.61 


630±155 


540.9 


62.4 


2.23±0.61 


564±155 


601.9 


59.3 


3.36±0.63 


850±160 


640.8 


18.4 


2.57±1.12 


651±284 


792.2 


19.5 


2.96±1.97 


749±277 


1005 


37.1 


4.93±0.83 


1248±209 


1044 


41.0 


2.48±0.75 


627±191 


1085 


41.1 


2.45±0.75 


620±190 


1126 


41.0 


2.48±0.75 


627±191 


5942 


200 


<0.99 


<253 


11729 


885 


<0.49 


<124 



B 


734.0 


10.0 


3.81±2.13 


309±173 




6044 


393 


<1.04 


<84 




18286 


615 


<0.83 


<67.5 



714.0 


9.8 


2.74±1.30 


168±80 


865.0 


9.8 


2.48±1.29 


152±79 


5839 


393 


<0.86 


<53 


17516 


885 


<0.56 


<34 



UVWl 


765.0 


39.8 


<1.47 


<143 


UVM2 


740.0 


39.8 


<0.84 


<104 


UVW2 


767.5 


19.0 


<1.68 


<139 



counterpart to the GRB, we performed an identical analysis 
on another star in the same field of view. This star reveals 
no change of flux throughout the observations, ruling out 
instrumental effects as the cause of the flux variations seen 
in the GRB OT. Observations carried out for a further 2 
days suggest the star superposed on the OT position can be 
considered a constant source, making it highly unlikely that 
the star in the XRT error circle flared at such time that it 
mimicked a GRB afterglow. 



6 SPECTRAL ENERGY DISTRIBUTIONS 

We constructed three optical to X-ray spectral energy dis- 
tributions (SEDs) spanning the prompt emission phase and 
one SED in t he afterglow phase using the ISIS spectral fit- 
ting package (JHouck fc Denicolall2000|) and fitting in coun t 
space using the method described in Istarling et al.l (|2007a ) . 
Throughout these fits we fixed the total X-ray absorp- 
tion (Galactic -(-intrinsic) to the value of A'h = 0.4 x 10^^ 
cm~^ derived from the X-ray spectral fits. Galactic ex- 
tinction is also included in all fits with E[B — V) = 0.4 
jSchlegel. Finkbeiner fc David 1 19981 ). We adopt the Small 
Magellanic Cloud (SMC) extinction curve for fits to the in- 



trinsic optical extinction in the host galaxjQ, allowing for 
distances to the source of 2: = 0-3 (see Section 5). 

The afterglow spectral energy distribution at 5670 s 
comprises UVOT B and V upper limits and XRT data. Fits 
to the SED are consistent with the synchrotron model, if 
there is a break in between the optical and X-ray bands 
(discussed further in Section 7.3). Fixing the difference be- 
tween the power law slopes in this broken power law model 
to 0.5, as expected for a cooling break, implies that some 
intrinsic absorption is also needed. Fitting to the optical 
upper limits and X-ray data we can set a lower limit to the 
required optical extinction of E{B — V)^ 0.1. This amount 
of intrinsic optical extinc tion lies at the high end of the 
observed distribution (e.g. lStarling et al.ll2007bl ). whilst the 
intrinsic X-ray absorpt ion is relatively small for a GRB (e.g. 
ICampana et al.ll2006ah . 

The optical emission observed at earlier times is not 
necessarily associated with the high energy prompt emission. 
In creating the SEDs for these prompt phases we assume the 
optical and X-ray/7-ray emission do come from the same 
component, and determine what spectral shape and amount 
of intrinsic extinction would be required in this scenario. The 
first of the prompt phase SEDs at 284 s includes UVOT V 
band, XRT and BAT data. We fit a Band function to the 
spectrum with the high energy photon index /3 fixed to 2.36. 
The availability of optical and high energy data provides 
a rare opportunity for Band function fits to a broadband 
SED. The V band datapoint lies below the best possible 
fit to the whole SED. To obtain an acceptable fit either 
intrinsic extinction would have to be present at the level of 
E{B — V) ^ 0.05 or both extinction and a further spectral 
break could be present. 

The second epoch SED at 770 s includes UVOT U and 

V band and XRT data. We flt an absorbed power law to 
the spectrum. The power law model overpredicts the optical 
data, illustrated in Fig. [S] and allowing for some intrinsic 
extinction decreases the Xreduced value by ~50%. In order 
to obtain an acceptable fit either intrinsic extinction would 
have to be present at the level of E{B - V) = 0.03-1.40, 
or we have a spectral break plus extinction. For the latter 
model we flnd a best fitting break energy of 0.7±0.2 keV, Fi 
= l-4ta2. r2 = 1.97±0.05 consistent with BAT-XRT Band 
function fits (Section 4.2) and E{B - V) ^ 0.5. The power 
law slopes differ by ~0.5 consistent with a cooling break 
inferring that the cooling break moves to lower energies with 
time. Both single and broken power law fits are, however, 
statistically acceptable and indistinguishable. 

The third epoch SED at 1082 s is the least constrained: 
this SED covers the flnal part of the steep X-ray decay phase 
of the light curve where the X-ray source has become fainter 
and XRT was operating in PC mode. A joint fit of the UVOT 

V band data and the X-ray spectrum manages to accommo- 



^ The SMC extinction curve is generally a better fit to GRB 
host galaxy extinction than that of the Milky Way or the Larg e 
Magellanic Cloud (e.g. IStarling et ahllJOOTbl : ISchadv et al.ll2007li . 
The SMC has the lowest metallicity of these three nearby galax- 
ies for which extinction curves have been derived, more closely 
resembling the low metallicity GRB host galaxies. We also note 
that Solar metallicity is assumed for the X-ray absorption mea- 
surements which means the adopted Nn intrinsic value is likely a 
lower limit. 
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Figure 8. The SED at 770 s after the trigger: XRT and V and U 
band data and UV filter upper limits. The unfolded (unabsorbed) 
model shown as a solid line is a power law fit to the observed V, 
U and X-ray data, showing the overprediction of the low energy 
spectrum. 
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Figure 9. The XRT light curve and fit to the energy injection 
time for the last occurring X-ray flare using the curvature model. 
A variable count rate to flux conversion was applied. This fit gives 



aTn 



of 632^ 



date the single optical point without inclusion of any intrin- 
sic extinction, with a slope of F = I.TGI^qq^ and x'^/dof = 
8.6/7, but is inconsistent with the fit to the X-ray data alone 
of r = 2.3 ± 0.3 (Table [2 Section 4.1) and with the after- 
glow SED fits. Including an intrinsic extinction component 
the power law photon index is then best fit with F = 2.3±0.3 
and the range of allowed extinction values is E{B — V) = 
0.5-2.2. More complex models cannot be tested with this 
dataset. 



7 DISCUSSION 

In GRB 070616 we have the rare opportunity to track the 
detailed spectral and temporal evolution from 7-rays to X- 
rays throughout the entire prompt emission phase. Hence we 
devote the majority of this discussion to the possible origins 
of the complex prompt emission, as well as the early optical 



emission, using the observational results given above. We 
also provide a brief discussion of the afterglow characteristics 
in the framework of the standard synchrotron model. 



7.1 The prompt emission mechanisms 

The prompt emission of GRB 070616 is atypical of GRBs 
in that the emission rises relatively slowly over about 100 
seconds to a peak, then persists at a fairly constant level 
for 150 s in 7-rays and 350 s in X-rays before showing a 
rapid decline. Throughout the ~ constant phase the light 
curve resembles a large number of fiares or fiickering su- 
perimposed on an underlying constant intensity emission. 
During this prompt phase we observe strong spectral evo- 
lution from hard to soft. The spectral evolution begins at 
To-f285 s at 7-ray energies, while the X-ray fiux is still at 
an approximately constant level, and begins 200 s later at 
X-ray energies around the onset of the steep X-ray decay 
(FigO. 

The spectral variability suggests a peak moving through 
the bandpass. We find that the peak of the spectrum can be 
accurately measured by adding the higher energy Suzaku 
WAM data to the Swift coverage. Assuming the Band func- 
tion to model the spectrum, the spectral peak started out 
at the high energy end of the BAT range — and was prob- 
ably above 200 keV before the XRT or Suzaku observations 
began — and moved to lower energies to most probably 
lie within the XRT bandpass beyond To+700 s (Fig. [5]). In- 
deed such strong spectral evolution may help explain some of 
the steep X-ray decline. We also find an indication that the 
BAT spectral lags may have decreased during the prompt 
emission from the first main peak to ~200 s later: such lag 
evolution appears to be common among spectrally evolving 
GRBs (J. Cannizzo et al. in preparation). The second, lower 
lag measurement was taken just after the BAT hardness ra- 
tio had begun its decrease and the spectral peak energy had 
fallen below 100 keV to be clearly detectable within the BAT 
energy range. 

While the spectral peak is moving we also measure a 
softening of the spectrum at frequencies below the peak. The 
spectral index above t he peak energy, 3, is f ound to be 2.36 
from the Suzaku data (jMorigami et al.ll2007l ). The spectrum 
is observed to be harder at higher count rates, as has been 
shown for several other GRBs (e.g. GRB 051117a, Goad et 
al. 2007a; GRB 061121, Page et al. 2007). Superimposed on 
this general trend is a hardening during the rise of individual 
flares, but this effect is not driving the overall evolution (Fig. 

Spectral evolution through the prompt phase has been 
noted previously, and is inconsistent with the idea that the 
curvature effect alone - the delayed arrival of emission from 
progressively higher latitudes within the jet - is driving 
the e mission during the steep decay phase (e.g. iLiang et al.l 
120061 ). If we try to force a curvature effect model to ex- 
plain the observed steep X-ray slope from the peak of the 
last flare at To-|-757 s out to To+950 s, we would expect 
that phase to begin at To+632tj2 s given the requirement 
of the model that a = 2 + f3 {F,^ oc t~°'v~^ and adopting 
the average spectral energy index over that time interval of 
P — 1.06, Fig. [9]). This would require a long initial emission 
period of over 600 s. In the sample of Swift GRBs studied by 
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Figure 10. Spectral evolution in the XRT band, plotted as f3 
vs. time for a direct comparison with Zhang et al. (2007), where 
/3 is the power law energy index of a single power law fit to the 
XRT spectra. We model this with the curvature model plus an 
underlying afterglow component to the strong spectral evolution 
(solid line). 



iLiang et al.l (|2006l ) the typical injection of energy tends to 
occur at the onset of the previous large peak. The curvature 
effect clearly cannot explain the steepest part of the X-ray 
decay from To-|-976-1200s (see last point in Fig.O also seen 
in e.g. GRB 051117a, Goad et al. 2007a; GRB 070110, Troja 
et al. 2007) . However, the combination of the curvature ef- 
fect and the strong spectral evolution we observe, which is 
a continuous process starting well before 600 s, may be able 
to a ccount for th e steep X-ray decay. 

IZhang et al.l l|2007l ') studied the X-ray tails of the 
prompt emission for a sample of Swift GRBs and found that 
only 25% of those could be fit with the curvature effect alone; 
those not well fit showed spectral evolution. From fits to 16 
GRBs, they tested and subsequently disfavoured two pos- 
sible causes of the spectrally evolving X-ray tails, namely 
an angle- dependent spectral index in structured jets and a 
superposition of the curvature effect and a power law de- 
cay component. The observed spectral softening could, they 
suggest, be caused by cooling of the plasma where the cool- 
ing frequency decreases with time. This manifests itself as 
a cut-off power law shape with the cut-off moving to lower 
energies with time, shown to be a goo d fit to GRB 06021 8 
and proposed earlier for GRB 980923 (JGiblin et all 1 19991 ). 
and similar to our Band function results for GRB 070616. 
In this scenario the peak energy we track through the BAT 
band would be the cooling frequency. The prompt emission 
SED at 770 s after the trigger, including optical emission, 
is consistent with an absorbed broken power law where the 
break shows the characteristics of a cooling break (though a 
single absorbed power law provides an equally good fit). In 
the afterglow SED we most probably do detect the cooling 
break, hence we may have observed it moving to lower ener- 
gies with time. However, we also observe softening of the low 
energy power law slope in addition to movement of a peak 
or cooling frequency. Oth er models for the prompt emission 
have been suggested - see IZhangj (|2007l ) for a review. 

We now investigate origins for the apparent spectral 
evolution in the simple superposition of multiple compo- 
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time since fitted peak (trigger+121 s) 

Figure 11. Top panel: UVOT V band flux light curve. Lower 
panel: BAT-XRT light curve with two-component model fit fol- 
lowing Willingale et al. (2007, solid line). This two-component 
model has proven a good fit to the majority of Swift high en- 
ergy light curves, but is unable to fit all the prompt emission of 
GRB 070616. The dashed line shows the flux proflle from the dou- 
ble power law model (equation 2) used to fit the /3 profile (Fig. 
IIOJ I. This model also cannot fit the observed light curve. 



nents. Fig. llOl shows the spectral evolution in the XRT band, 
characterised by the fitted power law slope /3 (/3 = F — 1), 
displa ying a similar beha viour to many of the bursts stud- 
ied bv lZhang et al.l (|2007|) and resembling most closely that 
of GRB 060510B which also had a lengthy prompt emission 
duration of T90 = 275 s with some fiaring (IBarthelmv et al.l 
I2OO61 I. To approximate the shape of the spectral evolution of 
GRB 070616 in the X-ray band (J3 vs. time) we fit a model 
with two power laws: 

-f^l _L 4_-f-"2 I7-/32 



FE = Ait-'''E-P'+A2t- 



'-E- 



(2) 



where t is time since a fitted To, At value and E is the pho- 
ton energy. We fixed the spectral slopes to those of the typ- 
ical measured prompt and afterglow values, /3i = 0.1 and 
/32 = 1.6. We fixed 02 at the probable afterglow temporal 
slope during the plateau of 0.5 (see Section 3) while cti , To, At 
and the relative contributions of the prompt and afterglow 
components were left to vary. This double power law model, 
overlaid on the evolution of /3 in Fig. 1101 appears to be a 
reasonable description of the behaviour through the steep 
decay and afterglow phases, with ai = 3.0, To^flt = To+200 
s and the two components contributing equally to the flux at 
To+800 s. The spectral shape evolves due to the increasing 
contribution from the afterglow to the overall spectrum. We 
then returned to the high energy light curve to test the dou- 
ble power law model fitted above. The dashed line in the 
lower panel of Fig. [TT] shows the double power law model 
applied to the high energy light curve. Clearly, although the 
evolution of (3 is well modelled the flux is not. We can there- 
fore rule out a spectrally invariant double power law model. 
It is possible that there is a further emission component, 
in addition to the high latitude emission, contributing to 
the steep decay ending the prompt emission, which becomes 
prominent between the prompt and early afterglow phases. 
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Table 4. A two-component fit to the BAT-XRT light curve fol- 
lowing Willingale ct al. (2007). The trigger time has been rescaled 
to To+121 s. Quantities Qp and Tp refer to the prompt compo- 
nent power law decay slope and start of that decay respectively, 
and O-a and Ta refer to the afterglow component, relative to the 
new zero time. 

component parameter value 



Ta 



7 07-1-0.12 
(.0(_Q 27 

14001?!^ s 

1.69«jg 

14380ti^« s 



To exa mine this possi b ility w e ap plied the modelli r iK pro ce- 
dure of lO'Brien et al.1 (|2006i ) and Iwillingale et all (|2007l ) in 
which it was shown that GRB early light curves can be well 
modelled by up to two emission components each consisting 
of an exponential followed by a power law decay. The first 
component models the prompt emission, adjusting the zero 
time. To, to provide the best fit. For GRB 070616 this moves 
zero close to the brightest peak in the prompt light curve 
at 121 seconds after the trigger. The best-fit model to the 
Swift BAT and XRT data is given in Table Hand Fig.fTT] 

We find that the light curve of GRB 070616 is not weU 
fitted by the usual two components consisting of an expo- 
nential plus a power law decline. The first, prompt compo- 
nent can fit the initial prompt data following the brightest 
peak, but then declines until the second, afterglow compo- 
nent starts to dominate. This final component can be fit with 
a power law of slope a ^ 1.7 from 4600 s onwards, typical of 
a post-plateau GRB afterglow (Table |4|. There remains an 
excess flux dominating from some 200-1000 seconds post- 
trigger. We speculate that this flux plus the spectral evolu- 
tion indicate an additional source of prompt emission which 
does not dominate in the majority of GRBs. 

Interestingly, several, although no t all, of the GRBs 
with strong spectral evolution studied in lZhang et al.l (|2007l l 
are also not well fitted by the two-compo nent light curve 
mode l (e.g. GRBs 051227 and 060614, see IWiUingale et al.1 
I2OO7I ). Others not well fitted by this mode lling procedure, 
su ch as GRB051117A IJGoad et all 1200731 ) were identified 
by IZhang et al.l (120071) as having evolution which they at- 
tributed to flares. Such GRBs may instead be more similar 
to GRB 070616. GRB051117A, for example, shows a long, 
bright, slowly declining prompt emission period with flaring 
followed by a sharp drop to a plateau or afterglow phase. 
GRB 070616 may be a case where the flux remains fairly con- 
stant, with flaring, up to the sharp drop. Among the Swift 
BAT bursts we notice of order 10 bursts with an initially 
flat high energy light curve, most prominent in GRB 070616, 
which when studied as a sample may provide more evidence 
in support of either a long-duration central engine or an 
additional component (N. Lyons et al. in preparation). 



T.2 The origin of the early optical emission 

The Swift era has provided a wealth of information on the 
early 7-ray and X-ray emission. The prompt optical emission 
has been harder to reveal as in most GRBs the flux level lies 



below that accessible to fast-response, modest aperture tele- 
scopes. Even so, for those GRBs which have been detected 
the relationship between the optical and X-ray emission is 
complex. In a very small number of cases the optical actually 
lies above the ext rapolation of the hi gher-energy emission 
(e.g. GRB 990123. lAkerlof et"aLlll999D . The optical can be 
dominated by a co mponent related to the prompt high en - 
ergy emission (e.g. IVestrand et al.ll2005l : lYost et al.ll2007al ). 
More usually, the optical lies below such an extrapolation, 
as shown for a sample of ROTSE-detected Swift GRBs by 
Yost et al. (2007b, using BAT data only for the extrapola- 
tions) . This pattern of diverse behaviour strongly suggests a 
mixture of components can contribute to the prompt optical 
emission, including the afterglow (due to the external shock) 
and emission from the prompt fireball (which may include 
a reverse shock). In the latter case the optical should con- 
nect, possibly with a spectral break, with the prompt 7-ray 
- X-ray spectrum. 

Few optical datasets have sufficient broadband data to 
check for consistency in spectral shape as well as overall 
flux level, limiting our ability to constrain reddening effects 
as opposed to spectral breaks. For GRB 070616 the combi- 
nation of long-duration prompt emission plus early optical 
multi-colour data can provide some constraints on the rela- 
tionship between the optical emission and the high-energy 
prompt emission. The V band light curve shows a rise to 
a plateau-like phase from the flrst observation beginning at 
To-f250 s and extending to To-f 1000 s. The optical decay 
does not follow the steeply decaying X- and 7-ray prompt 
emission. A gradual rise would be typical of the start of af- 
ter glow emission, seen in a number of bursts (jMolinari et al.l 
[20071). We performed a power law fit to the V band light 
curve shown in Fig. [T]and measure a = -(0.7 ±0.3). This is 
consistent with predictions for the afterglow rise in the fire- 
baU model of a = -0.5 (e.g. Zhang & Mcszaros 2004). Fits 
to the three prompt emission SEDs for GRB 070616 all show 
that extrapolation of the high energy spectrum to lower en- 
ergies overpredicts the optical flux, however this may be 
overcome with a relatively large amount of intrisic extinction 
or a combination of extinction and a spectral break, if the 
optical data are assumed to be prompt emission dominated. 

We can compare this burst with other GRBs for which 
well-time-resolved prompt optical data are available, and 
whose light curves are not afterglow-dominated from early 
times. For GRB 050820A iVestrand et al. (2006 ) find a ris- 
ing optical component plus a component correlated with 
the 7-ray emission. In GRB 050904 the optical data rose 
to a plateau before mimicking a fiare seen in the X- 
ravs (iBoer et ahlbOOd ). Similarly in GRB 051111 lYost et all 
(|2007al ) suggest a combination of components produce the 
prompt optical emission with a significant contribution to 
the early optical fiux from an extension of the higher-energy 
spectrum. These are similar to the case of GRB 061121 
(|Page et al.l 120071 ) where the high energy prompt emission 
consisted of a single large peak, which was mimicked in the 
optical data, although an additional spectrally evolving op- 
tical component is also seen. Extrapolation of the high en- 
ergy spectrum in GRB 061121 overpredicts the early optical 
emission, suggesting a spectral break. We note that the ex- 
trapolation of the high energy spectrum to the optical band 
predicts a relatively higher optical flux in GRB 061121 than 
in GRB 070616, due to its softer high energy spectrum. 
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For GRB 070616 we conclude that strong optical emis- 
sion associated with the prompt phase and described by 
a simple extrapolation of the observed higher energy spec- 
trum cannot be present without invoking a further spectral 
break in between optical and X-rays or extinction in the 
host at the highest end of the currently o bserved values (e.g. 
ISchadv et ai]|2007l : IStarling et al.ll2007bl ). We can likely rule 
out t he "external-external" shock model (iMe szaros fc ReesI 
I1993I ). which requires that the broadband prompt emission 
is produced by one and the same mechanism and any vari- 
ability is caused by density enhancements in the external 
medium, given the very different behaviour of the optical 
and X-ray light curves. 

7.3 The afterglow 

After To +4000 seconds the X-ray afterglow can be consid- 
ered 'typical' of what we know for GRB afterglows: the 
spectral and temporal parameters lie within the range of 
the ~250 Stmft-ohserYeA GRBs. We test these parameters 
against the synchrotron model, using the closure relations 
for the case of a cooling break in between the optical and 
X-ray bands, Vc < vx and a = (SF — 4)/2, and find that the 
observed temporal and spectral X-ray slopes are just consis- 
tent at the 90% level and therefore follow the expectations 
for synchrotron emission (in the slow cooling regime). This 
result also indicates that the afterglow has not reached the 
'jet' phase, in agreement with the lack of evidence for a jet 
break during the Stmft observations. We do not find good 
agreement if the cooling frequency were to lie above the 
X-ray band, Uc > vx, at this time, consistent with the re- 
quirement of a broken power law to best describe the SED 
of the late-time afterglow spectrum plus optical upper lim- 
its. Therefore we conclude that the observed spectral break 
is due to the cooling frequency. The cooling frequency is 
hkely to lie at or above the X-ray frequencies at earlier times 
(<4000 s), depending on the circumburst medium density 
structure, but we cannot confirm this since the initial X-ray 
emission is prompt-dominated. 

We find a 90% confidence range for the injected elec- 
tron energy index, p, of 2.68-3.46, as determined from the 
measured X-ray spectral index where F = (p + 2)/2. This 
is higher than the commonly quoted 'universal' value of 
p = 2.2 derived from some numerical si mulations of parti- 
cle ac celeration in relativistic shocks (e.g. lAchterberg et al.l 
I2OOII ). but other numerical simulations (e.g. iBarind |2004| ) 
and recentl y publi shed observational studies of sample s of 
afterg lows IjShen. K umar & Robinson 2006; Starling c t al.l 
[2007g) find a range of allowed values for p which encom- 



pass the values we report for GRB 070616. 



8 CONCLUSIONS 

GRB 070616 has afforded us the rare opportunity to track 
the detailed spectral and temporal evolution from 7-ray to 
optical wavelengths throughout much of the prompt emis- 
sion phase, owing to Swift's fast slew capability, the broad- 
band coverage provided by the combination of Swift with 
Suzaku WAM, and the long prompt emission duration. The 
high energy light curve remains generally flat for several 
hundred seconds before going into a steep decline. 



Spectral evolution from hard to soft is clearly tak- 
ing place through the prompt emission (flat and decaying 
phases) of GRB 070616, beginning at 285 s out to at least 
1000 s. The high energy spectrum from 0.3-800 keV is well 
modelled by a Band function. We track the spectral peak 
energy moving from 135 keV down to ~4 keV and measure 
a softening of the low energy spectral slope from F~l.l- 
2.3. The curvature effect, whilst contributing to the spectral 
shape, is not dominating during the steep decay at ~1000 s. 

The presence of an additional component, perhaps not 
present in most GRBs, is favoured by the inability of a two- 
component fit to model the light curve. This unusual light 
curve shape of an initial constant flux followed by a sharp 
drop, plus the strong spectral evolution, is most prominent 
in GRB 070616, but we note a similar structure can be found 
in a small number of other GRBs which are also not dom- 
inated by the curvature effect and not well fit by the two- 
component model. 

We find that the optical data during the prompt phase 
are not consistent with a simple extrapolation of the high 
energy spectrum, requiring either significant intrinsic extinc- 
tion or some extinction plus a break in the spectrum between 
the X-ray and optical bands. 

The afterglow is consistent with the synchrotron model, 
with the cooling break lying in between optical and X-ray 
bands at a few thousand seconds after the trigger. 

The prompt emission of GRB 070616 comprises a com- 
ponent well fitted with a Band function and a possible fur- 
ther component. It is clear that both broadband coverage 
and good time resolution are crucial to pinning down the 
origins of the complex prompt emission in GRBs. 
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Table 2 Joint spectral fits, using the energy ranges XRT 0.3-10 keV; BAT 15-130 keV and WAM 150-800 keV, and over varying time ranges covering 137—963 s since 
trigger. Three models are tested; power law, broken power law and Band function. All spectral slopes reported here are photon indices. 
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